First-Principles Calculation of the Superconducting Transition in MgB2 within the 

Anisotropic Eliashberg Formalism 
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We present a study of the superconducting transition in MgB2 using the ab-initio pseudopotential 
density functional method and the fully anisotropic Eliashberg equation. Our study shows that 
the anisotropic Eliashberg equation, constructed with ab-initio calculated momentum-dependent 
electron-phonon interaction and anharmonic phonon frequencies, yields an average electron-phonon 
coupling constant A — 0.61, a transition temperature Tc = 39K, and a boron isotope-effect exponent 
Qfl = 0.31 with a reasonable assumption of fi* = 0.12. The calculated values for Tc, A, and a_g are in 
excellent agreement with transport, specific heat, and isotope effect measurements respectively. The 
individual values of the electron-phonon coupling A(fc, k') on the various pieces of the Fermi surface 
however vary from 0.1 to 2.5. The observed Tc is a result of both the raising effect of anisotropy 
in the electron-phonon couplings and the lowering effect of anharmonicity in the relevant phonon 
modes. 



Although MgB2 is a readily available sp-bonded mate- 
rial, superconductivity in this material with a transition 
temperature of Tc = 39 K was found only very recently 
0. This relatively high Tc has motivated many stud- 
ies, as has the observation that the detailed supercon- 
ducting properties of MgB2 show significant deviations 
from those calculated using the standard BCS model. 
The isotope effect exponent for boron is reduced sub- 
stantially from the conventional value for sp metals , 
and the average electron-phonon coupling strength A ob- 
tained from specific heat measurement |^,^ seems too 
small to justify the high Tc. In addition, specific heat 
measurements tunneling [|j and photoemission |^] 

spectra, and point-contact spectroscopy ||J^ show low 
energy excitations suggesting a secondary gap. Theoret- 
ical calculations show that the Fermi surface has sev- 
eral pieces and is very anisotropic ||l^, and that the 
electron-phonon coupling is dominated by the in-plane 
B-B stretching modes (i?2g) ||lC|-[l^ which have a large 
anharmonicity The electron-phonon interaction 

varies strongly on the Fermi surface jlj,|l^], and a two- 
band model suggests a multigap scenario |l4||- However, 
there has not yet been a quantitative, first-principles cal- 
culation of Tc including the full variation of the electron- 
phonon interaction on the Fermi surface and the anhar- 
monicity of the phonons to help confirm the phonon- 
mediating pairing mechanism for superconductivity in 
MgB2. 

In this letter, we present Tc and isotope-effect expo- 
nents for MgB2 obtained by solving the k and uj depen- 
dent Eliashberg equation. It is shown that the anisotropy 
(i.e., the electronic-state dependence) of the electron- 
phonon interaction on the Fermi surface is strong enough 
to raise to 39K even though the interaction is weak- 
ened by the anharmonicity of the phonons as compared 
to the harmonic case. In addition, it is shown that the 
anharmonicity of the phonons reduces as to 0.31. These 
results show that conventional phonon- mediated electron 



pairing theory can explain superconductivity in MgB2 
when both the anisotropy of the electron-phonon inter- 
action and the anharmonicity of the phonons are properly 
taken into account. The solution of the full Eliashberg 
equation at low T further yields different gap values for 
the different parts of the Fermi surface. The gap value 
distribution clusters into two groups - a small value of ~ 
2 meV and a large value of ^ 7 meV. This feature and its 
physical consequences will be described in more details 
in a future publication . 

The phonon frequencies and electron-phonon matrix 
elements are calculated using ab-initio pseudopotentials 
and the local density approximation. We used a 18 x 
18 X 12 k-point grid in the Brillouin zone (BZ), and in- 
cluded planewaves up to 60 Ry as a basis to expand the 
electronic wavefunctions. The calculated equilibrium lat- 
tice constants are a = 3.071A and c = 3.578A, m good 
agreement with measured values [|l| . We performed total 
energy calculations with frozen phonons for all phonon 
modes at all the high symmetry points of the BZ. The 
variation of the total energy with a frozen phonon am- 
plitude is fitted with a fourth order polynomial to ac- 
count for the phonon anharmonicity. To obtain the har- 
monic phonon frequency, we use the quadratic term of 
the fitted curve and calculate the frequency classically; 
whereas for the anharmonic phonon frequency, we cal- 
culate quantum-mechanical vibrational states including 
the anharmonic terms, and take the energy difference 
of the two lowest states. In the case of the degenerate, 
in-plane B-B stretching modes (i?2g) at F and A, we 
calculate quantum-mechanical vibrational states in two 
dimensions after the total energy is fitted in a plane with 
E{r,9) = Eo + C2r^ + Cir^ + {c3r^-\-C5r^)cos{39). We use 
natural atomic weights for B and Mg, that is, 10.81 for 
B and 24.31 for Mg, but ^"B or ^^Mg are used when we 
recalculate the phonon frequency for the isotope effect. 
The linear electron-phonon matrix elements are evalu- 
ated directly from the total self-consistent change in the 
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crystal potential caused by a frozen phonon. 

Table I shows the frequency of the in-plane B-B 
stretching mode (£'29) at T. This mode is doubly degen- 
erate along the line from F to A, and has a large anhar- 
monicity and a large electron-phonon coupling. Anhar- 
monicity increases the frequency by 20% and weakens the 
corresponding electron-phonon couplings by 30%. The 
calculated anharmonic frequency, 75.9 meV, for the i?2g 
mode at F agrees very well with the results from Raman 
measurements (75.9 meV and 76.9 meV [|8|) as well 
as other theoretical calculations The E2g modes 

at M, L, K, and H however have very little anharmonic- 
ity and small electron-phonon coupling. The strong an- 
harmonicity and the large electron-phonon coupling are 
thus confined to phonons in a small volume in k-space 
near the F to ^ line. 

The calculated phonon frequencies and electron- 
phonon matrix elements — {k\SVi\k') for the jth 
phonon mode are interpolated onto a 18 x 18 x 12 grid in 
the BZ through the following three-step process. First, 
we interpolate the dynamical matrices using a weighted 
average of those at the symmetry points and obtain the 
phonon frequencies and eigenvectors on the grid by di- 
agonalizing the dynamical matrices. Second, we inter- 
polate the induced crystal potential change by a phonon 
on the grid from the calculated crystal potential changes 
at the symmetry points using weighting factors deter- 
mined from the phonon frequencies and polarization vec- 
tors calculated on the fine grid. Finally, we calculate 
the electron-phonon matrix elements on the grid using 
the interpolated crystal potential change. All calcula- 
tions are done twice for comparison: one with harmonic 
phonon frequencies and another with anharmonic phonon 
frequencies. To study the isotope effect, we repeat the 
entire procedure with an isotopic atomic mass. 
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FIG. 1. Phonon density of states F{uj) and the isotropic 
Eliashberg function a'^F{uj) for MgB2. 

Figure 1 shows the phonon density of states F{uj) 
and the standard Eliashberg function a^Floj). The 
phonon density of states shows a large peak at 37 meV 
arising from the van Hove singularities in the acoustic 
phonons, but these phonons make no significant contri- 
bution to a^F^uj). There is a large dominant peak in 
a'^F{LLj) at 63 meV for the case of harmonic phonons 
but at 77 meV for anharmonic phonons. The dom- 
inant peak in a^F{uj) is caused by the in-plane B-B 



stretching modes {E2g). Because the i?2g modes are 
highly anharmonic and have very large electron-phonon 
coupling only for phonons within a small volume along 
the F to A line in k-space, anharmonicity has little ef- 
fect on but it causes a big shift in a'^F[ijj). In 
the case of harmonic phonons, as is shown in Table I, 
the isotropic average electron-phonon coupling constant, 
A = 2 / du!a'^F{uj)/uj, is 0.73 and the logarithmic average 
frequency, cj;„ = cxp[(2/A) / da;a^F(u;) Incj/w], is 59.4 
meV. These values and the overall shape of a^F{uj) with- 
out anharmonicity in the present calculation are in good 
agreement with previous calculations With an- 
harmonicity, A is reduced to 0.61 and ujin is increased to 
63.5 meV. Since A corresponds to the mass enhancement 
factor for the density of states at the Fermi level regard- 
less of anisotropy in the electron-phonon interaction , 
we can compare the calculated A with results of specific 
heat measurements. The reduced value of A = 0.61 due 
to anharmonicity agrees very well with result of specific 
heat measurements which give a A of 0.58 § and 0.62 [|. 
This agreement is evidence that phonon anharmonicity 
weakens the electron-phonon interaction in MgB2 . How- 
ever if this value of A = 0.61 is used in the McMillan |2^] 
or the Allen-Dynes formula for Tc, the predicted Tc 
would be far lower than experiment. 

Unlike previous studies, we solve the fully anisotropic 
Eliashberg equation for superconductivity in MgB2. The 
anisotropic Eliashberg equation at Tc [M is 
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and Wr is 

A; 



X A(fc',iw„') , 

where w„ = {2n + 1)t:Tc, /„ = l/\2n + 1 
the fraction of the density of states at k on the Fermi 
surface. For the definition of Z, A, X{k,k' ,n), and s„, 
see Ref. With the exception of our calculation 

of the phonon frequencies and electron-phonon interac- 
tion provides all the material parameters for solving Eq. 
(1) and hence for obtaining Tc from first principles. The 
dimensionless Coulomb pseudopotential, /x* ||2^, is the 
only free parameter in our calculation; but it is known 
to be of order 0.1 in most metals 22 2^ and we show 
below that the superconducting properties of MgB2 are 
insensitive to fi* . For comparison, we also calculate Tc 
using the isotropic Eliashberg equation. 
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where X{n) = J2kk' Wj^Wj:, X{k , k' , n) . Hence A(n) is the 
electron-phonon coupling averaged over all pairs of {k, k') 
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on the Fermi surface. (A(n = 0) is equal to the spe- 
cific heat A discussed above.) The isotropic Eliashberg 
equation is thus a special limited case of the more gen- 
eral anisotropic equation. If the electron-phonon interac- 
tion A(fc, fc', n) did not depend strongly on the electronic 
states on the Fermi surface, the isotropic equation would 
be an appropriate approximation. 




FIG. 2. (color) Variation of the electron-phonon interac- 
tion A on the Fermi surface of MgB2 . The left plot shows the 
mass enhancement factor given by A(fc, n — 0). The right plot 
shows A(fc = fco, fc', n = 0) as a function of k' for a fixed fco on 
the Fermi surface near F. 
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FIG. 3. (color) Number density of (fc, k') pairs on the Fermi 
surface versus the value of X{k,k' ,n = 0). The number den- 
sity is split into three sets: both k and k' are on the two a 
cylindrical sheets of the Fermi surface (red line) , both on the 
two TT tubular sheets (green line), and one on a a cylindrical 
sheet and the other on a tt tubular sheet (blue line). 

Figure 2 shows the variation of the calculated electron- 
phonon interaction on the Fermi surface of MgB2. The 
Fermi surface of A/IgB2 consists of four sheets; two hole- 
like coaxial cylinders consisting of cr-boron in-plane states 
along F to A, a holelike tubular network of 7r-boron 
states connecting regions near K and Af , and an electron- 
like tubular network of 7r-boron states connecting regions 
near H and L. In Fig. 2, the mass enhancement factor for 
states at k given by X{k,n = 0) — J2k' W^,X{k,k' ,n — 0) 
shows two well-separated sets of values on the Fermi sur- 
face. X{k, n — 0) is about 0.8 ~ 1.0 on the two a cylindri- 
cal sheets and is only about 0.3 ~ 0.5 on the two tt sheets. 
For more detail, we depict the value of A(fc = ko, k' , n = 0) 
as a function of k' for a fixed fc = fco on the Fermi surface 
near F. It shows strong and varying strength for scat- 



tering onto the a cylindrical surfaces but rather weak 
strength for scattering to the tt sheets. Figure 3 shows 
the number density of (fc, fc') pairs on the Fermi surface 
plotted as a function of the value of A(fc, fc', n = 0). The 
coupling strength A(fc, fc',n = 0) between states on the a 
cylindrical sheets has values exceeding 2.0 which is much 
larger than those within the tt tubular sheets or between 
a a cylindrical sheet and a tt tubular sheet. All this in- 
formation and A(fc, fc',n) with non-zero n are taken into 
account in solving the anisotropic Eliashberg equation. 

The anisotropic Eliashberg equation including anhar- 
monicity in the phonon frequencies yields 42K > Tc > 
37K for 0.10 < ^J* < 0.14. In particular, Tc is 39K when 
yU* = 0.12. To investigate the role of anisotropy in the 
electron-phonon interaction and of anharmonicity of the 
phonons, we calculate Tc disregarding one or the other, 
as shown in Table I. If we neglect the anisotropy and 
calculate with the isotropic Eliashberg equation, 
drops to 19K for — 0.12. This shows that the strong 
variation in the electron-phonon coupling of scattering 
on the Fermi surface is crucial to the observed high Tc in 
MgB2. As another comparison, if we calculate Tc using 
the anisotropic Eliashberg equation but neglecting the 
anharmonic effect in the phonon frequencies, Tc goes up 
to 55K for fi* = 0.12. Hence anharmonicity lowers Tc 
in MgB2. Thus we conclude that anisotropy in MgB2 is 
essential to produce the anomalously high Tc, especially 
in view of the fact that the electron-phonon interaction 
is weakened by anharmonicity. We note that in MgB2, 
an average electron-phonon coupling A cannot be cor- 
rectly determined from Tc using the McMillan (2^j2j] or 
Allen-Dynes equations |^l|, however a determination of 
A from the specific heat measurement is still valid. The 
Tc of MgB2 is not a function of the usual isotropically 
averaged electron-phonon interaction A given above; it 
depends on the details of electron-phonon interactions 
on the full Fermi surface. This explains the apparent 
discrepancy between the values of A estimated from spe- 
cific heat measurements and A estimated from Tc using 
simplied isotropic models. 

To calculate the isotope-effect exponent a {Tc cx 
M~"), we recalculate Tc using the mass of either ^°B or 
^^Mg in place of the natural atomic weight. Since the rel- 
evant phonons are related to motion of the B atoms, the 
Coulomb pseudopotential /i* also depends on the atomic 
mass of boron Mb- The change of /i* is simply 

^.• = -(.-)^gg, (3) 

using /i* — fi/[l + iJ,\n{ep/uJo)] and uja oc A/^°'^. 

Table I shows calculated isotope-effect exponents with 
fi* = 0.12 both with and without phonon anharmonicity. 
Without anharmonicity, the slight deviation of the sum 
of the two exponents, as and aMg, from the value of 
1/2 is due to the change of ^* given by Eq. (3). If we 
neglect Sfi*, the total isotope-effect exponent would be 
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0.5. In contrast, when we include anharmonicity in the 
phonon frequency, the isotope-effect exponent for boron 
is substantially suppressed. We obtain as — 0.31 and 
ctMg = 0.05 from the anisotropic Eliashberg equation 
with anharmonic phonon frequencies. In this case, the 
contribution of (5/x* to the decrease of as is only 0.02, so 
the anomalously low isotope-effect exponent is primarily 
due to phonon anharmonicity. 

In conclusion, we have shown from first-principles cal- 
culations that MgB2 is a conventional phonon-mediated 
superconductor whose properties require, for a correct 
description, a solution of the fully anisotropic Eliash- 
berg equation including phonon anharmonicity. The 
isotropic Eliashberg equation seriously underestimates Tc 
because it fails to account for the (fc, A:')-dependency of 
the electron-phonon interaction on the Fermi surface . 
We show that the electron-phonon coupling is exceed- 
ingly strong for certain pairs of (k, k') on the discon- 
nected Fermi surface of this material. The anisotropy of 
the electron-phonon interaction in MgB2 is strong enough 
to produce the observed Tc of 39 K in spite of a moder- 
ate average electron-phonon interaction as also seen in 
specific heat measurements. In addition, we have shown 
that the anharmonicity of the phonons in MgB2 weakens 
the electron-phonon interaction and reduces the boron 
isotope-effect exponent. 
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TABLE I. Transition temperature Tc and isotope-effect ex- 
ponents a with /i* = 0.12. Numbers in parentheses are the 
values of Qs when of Eq.(3) is ignored. The averaged elec- 
tron-phonon coupling A and the frequency coph of the in-plane 
B-B stretching modes (-E29) at P are also included. 
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